Spin Driven Jahn- Teller Distortion in a Pyrochlore system 
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The ground-state properties of the spin-1 antiferromagnetic Heisenberg model on the corner- 
sharing tetrahedra, pyrochlore lattice, is investigated. By breaking up each spin into a pair of 1/2- 
spins, the problem is reduced to the equivalent one of the spin- 1/2 tetrahedral network in analogy 
with the valence bond solid state in one dimension. The twofold degeneracy of the spin-singlets 
of a tetrahedron is lifted by a Jahn- Teller mechanism, leading to a cubic to tetragonal structural 
transition. It is proposed that the present mechanism is responsible for the phase transition observed 
in the spin-1 spinel compounds ZnV^CU and MgV204. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.80.+q 



Geometrically frustrated spin systems have been a fas- 
cinating subject to study since the Anderson's pioneering 
work on classical spins with a disordered ground state |Q . 
The pyrochlore and the spinel compounds, both including 
the three-dimensional (3D) tetrahedral network sharing 
the vertexes, are the typical examples of such systems in 
nature. As recent theoretical progress on pyrochlore spin 
systems, it may be mentioned that the spin- 1/2 antifer- 
romagnetic (AF) Heisenberg model is shown to have a 
spin liquid ground state [Q] and that, even for ferromag- 
netic Ising model, nontrivial effect of frustration leads to 
a newly found spin ice ground state ||. In this letter, we 
investigate the ground-state properties of the spin-1 AF 
Heisenberg model on the pyrochlore lattice. 

We follow a general strategy of constructing effective 
variational wave functions, which approximately describe 
the low energy states in the same spirit as the resonat- 
ing valence bond (RVB) Q or the one-dimensional (ID) 
valence bond solid (VBS) approaches ||. These simple 
pictures are of great use to understand essential physics of 
the systems in more intuitive way. Our simple scenario, 
where the low energy effective Hilbert space is assumed 
to adiabatically continue to the manifold of the product 
wave functions of the twofold degenerate spin-singlets at 
each tetrahedron, leads to a spontaneous breakdown of 
the lattice symmetry. This new effect of the magneto- 
elastic interaction in the pyrochlore spin system gives a 
consistent picture with the recent experimental results 
on the insulating spin-1 spinel compounds ZnV2 04 and 
MgV 2 4 

First, let us introduce the AF Heisenberg model on the 
pyrochlore lattice, Fig. |l|, in two different but equivalent 
representations where the interactions are defined for the 
sets of bonds and tetrahedra, respectively, 



const., (J > 0), 



<1,3> 



k=l 



where S is the spin-1 operator and < i,j > denotes a 
nearest neighbor pair. The index k, numbered from 1 
to AT 4 , specifies a tetrahedron and St ot represents the 



sum of the four spins on the k-th tetrahedron. For the 
present spinel compounds, we expect that the relatively 
small spin value makes the isotropic limit a good starting 
point. On the other hand, it is well known that there are 
strong anisotropies for some rare earth titanates [^) . 




FIG. 1. The ideal locations of the vanadiums and oxygens 
in ZnV204 and the network of the vanadium ions viewed from 
the [111] axis. 



To tackle this problem, we develop an analogous 
method with that of the VBS state, which captures essen- 
tial physics of the Haldane phase in ID spin-1 systems. 
By breaking up the original spin-1 into two spin- 1/2 ob- 
jects (denoted by s), we rewrite the Hamiltonian neglect- 
ing a constant energy shift as follows, 



JYi 



H = V 



■v, 



(2) 



fe=i 



where (s^, Sfe 2 , Sfc 3 , Sfc 4 ) are the spin-1/2 operators form- 
ing a tetrahedron. V is the operator to symmetrize the 
states spanned for all pairs of the 1/2-spins on every ver- 
tex and thus restricts the expanded Hilbert space to that 
(l)of the original spin-1 model. The four times larger cou- 
pling constants are required so as to describe the original 
exchange interactions by those within the 1/2-spins on 
each tetrahedron. The tetrahedron-unit representation of 
the spin-1 pyrochlore system has a significant difference 
from the ID VBS state in the sense that our fundamental 
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unit is not a bond with two 1/2-spins but a tetrahedron 
with four 1/2-spins jl2). In particular, there are twofold 
ground-state degeneracy which makes the situation more 
interesting as discussed below. We believe that this sim- 
ple transformation is useful to grasp the essence of the 
low energy physics of the original problem. 

Following the general recipe for the ID VBS state, we 
first discuss properties of the fundamental unit. Since the 
total spin of the tetrahedron is a good quantum number, 
the two spin-singlets constitute the degenerate ground 
states. By using the linearly dependent three singlets 
whose total sum equals zero, |2) = 1 1234) , |3) = |1342), 
and 1 4) = |l223*), where \ahcd) is the product of the two 
valence bonds, \ab) = (talfe — latfe) /V% an( l similarly 
defined \cd), we express the tetrahedron singlets in or- 
thonormal bases, 

|i?} = ^(|2)+c|3}+c*|4)), and \L) = \R)*, (3) 

where uj = exp(27ri/3) and oj* its complex conjugate. 
These are the bases of the E representation of the T<± 
group and are the eigen states of the ±27r/3 rotation 
around the four different trigonal axes. 

As a next step, we generate the direct products of the 
local singlets. By symmetrizing the two 1/2-spins on ev- 
ery vertexes, we obtain the states 

Ni 

\^ gs )=V]l(r k \Rk)+lk\L k }), (4) 
fe=i 

defined for arbitrary {rk,h} with \rk\ 2 + \lk\ 2 = 1. These 
states are the ground states of the Hamiltonian repre- 
sented by a sum of the projection operators onto the total 
spin three (represented by P 3 ) subspace for the tetrahe- 
dron made of four 1-spins. Because the four 1/2-spins, 
out of the broken-up eight 1/2-spins, form a singlet state, 
then the maximum total spin of the original tetrahedron 
must be less than three. To be explicit, the projection 
operator P 3 of the fc-th tetrahedron is given by, 

*-h < s: ") 2 -tSo ( s "") 4 +t ^ W-^s (*•)' m 

We assume that the Hamiltonian (|l|) and J^k P^ be- 
long to the same universality class and share the essen- 
tial properties in the same way as the relation between 
the AF Heisenberg model and the VBS model in one 
dimension JD|. This assumption leads to a consistent 
picture for the properties of ZnV2C>4 as minutely dis- 
cussed in the latter part of this letter. In the above 
argument, we supposed that the symmetrized states do 
not vanish. As a matter of fact, one can show by us- 
ing the Schwinger boson representation that the direct 
products of any linear combination of the form, eq.(Q), 
have the same nonzero norm even after the symmetriza- 
tion. Accordingly, although the exact dimension of the 
symmetrized states is not known, a macroscopic ground- 
state degeneracy of order ~ 2 Ni is expected for the model 



given by the projection operator. It is the key assump- 
tion of the present study that the low energy part of the 
Hilbert space of eq.(|l]) adiabatically continues to these 
symmetrized states of the singlets' products. 

In real materials, the ground-state entropy must be 
zero and it is quite natural to expect that the above- 
mentioned macroscopic degeneracy is lifted. In what 
follows, therefore, we concentrate on finding some rea- 
sonable mechanism which stabilizes the observed physi- 
cal state among a number of nearly degenerate ground 
states. As long as the lattice symmetry remains to be 
cubic, it would be appropriate to treat the problem as a 
pure spin system and investigate the effects of the longer- 
range interactions as well as the lattice topology as a 
source of lifting the degeneracy. One can argue that 
this type of lifting produces a kind of chiral ordering, 
which will be published elsewhere jyfl . When the inter- 
action between the lattice and spin degrees of freedom is 
more important, another way of lifting would take place 
with spontaneous breaking of the cubic lattice symmetry 
p5| . The energy loss from the lattice rigidity and gain 
from the magneto-elastic interaction determine the sta- 
ble structure, a generalized Jahn- Teller effect of a novel 
type. Thus it has some similarity with the spin-Peierls 
transition usually discussed in ID. 

For this purpose we consider a single tetrahedron with 
four vanadium sites (V1-V4) as shown in Fig. [|. To 
describe lattice vibrations around the stationary points 
(Ri), let us define small deviations by R4 — R- + Axi. 
According to the symmetry of the tetrahedron (Td point 
group), the normal modes are classified into the A\ (Qa- 
its normal coordinate), E (Q u ,Qv), and T 2 (Q^Q^jQc) 
representations after eliminating the uniform translation 
(T2) and the uniform rotation (Ti). The normal coordi- 
nates of Ai and E representations are written as, 

Q A {Ai) = (V + V + Z)/V3, (6) 
Q U (E) = (X- Y)/V2, (7) 
Q V (E) = (X + Y -2Z)/V6, (8) 

where X, Y, and Z are the uniform elongation for x,y, 
and z directions, given by (—Ax\ + Axi + Ax-& — Ax4,)/2, 
(-Ay! + Ay 2 - Ay 3 + Ay 4 )/2 and (-A Zl - Az 2 + Az 3 + 
Azi)/2, respectively, see Fig.|[ 




Qa(Ai) Q u (E) Q v (E) 

FIG. 2. The schematic representations of the normal modes 
viewed from the c axis. 
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From simple group theoretical consideration, the local 
Hamiltonian, Hk, for a single tetrahedron is given, up to 
the first order of deviations, by 



H k = V6J s a + g a 
+ 9t 



dJ 



OQa 
dJ 



QaSa+ 9et^t- 
"WE 



6Qt 3 



(Q 5 s ? + Q v s v + Qc s c) 



(9) 



where sr's are the bases of the irreducible representa- 
tions made from the bilinear combinations of the spin 
operators, 
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(10) 



To illustrate the coupling constants, let us assume for 
simplicity that the exchange coupling depends only on 
the distance between the spins. Then J — J(r) and 
the derivatives by each normal coordinate are given by 
2dJ/dr (Ai mode), dJ/dr (E), and VldJ/dr (T 2 ), re- 
spectively. 

In the ground-state spin-singlet subspace, the matrix 
elements for the T 2 modes vanish and thus may be ne- 
glected. By including the lattice restoring energy, Hk is 
given by 



V3gE dJ 



2 dQ E 



(Q u a x + Q v a z )+ ^ (Ql+ Q 2 ) 



(11) 



Here we have neglected the Qa mode which simply 
renormalizes the Heisenberg coupling constant to J + 
(dJ/dQA)QA/VQ- We have used the real basis, \u) — 
(\L) + \R)) /\/2 and \v) = (\L) - \R)) /V2i and the Pauli 
matrices are defined by a x = \u)(v\ + |i>)(it| and <r 2 = 
— Concerning the symmetrization V, its 

effect may be taken as a renormalization factor of the 
order of unity (gA, gE, and gx in eq(^)). In other words, 
the form of the local Hamiltonian is determined by the 
symmetry and the important properties such as the equal 
coupling constants for the two components of the E mode 
and the vanishing of the matrix elements of the s^, s^, 
and are generic independent of the symmetrization on 
every vertexes. 

Up to this order, there is no force to determine 9 of 
(Q v ,Qu)=(pcos9,pshi8), withp = V3g E \-^\/2k 2 E . In 
order to fix this phase 9, it is necessary to take some 
anisotropies into consideration, such as the anharmonic- 
ity of the lattice restoration energy |l6| or the higher or- 
der deviations of the magnetic exchange interaction p7[ . 
After some calculations including up to the third (for the 
lattice part) and fourth order (the spin part) of the lattice 



deviations, respectively, we find both terms to be propor- 
tional to cos 39, which means that the elongated or com- 
pressed lattice structure along one of the a, b, or c-axis 
becomes the lowest states depending on the sign before 
the cos 36* term. We must distort each tetrahedron coop- 
eratively all over the crystal without contradictions. For 
example, a uniform local tetrahedral distortion causes a 
uniform lattice metamorphosis and an alternating lattice 
distortion develops another type of lattice distortion. In 
any way, the many-fold degeneracy is lifted at this struc- 
tural transition and a quenched local tetrahedral struc- 
ture is realized. 

Now let us discuss an application of the present sce- 
nario to the spin-1 normal fee spinel compounds ZnV204 
[g-|l0| or isostructural MgV 2 4 where the mag- 

netic V 3+ ion, with (3d) 2 electronic state forming 5=1, 
resides on every lattice point of the pyrochlore lattice. 
These two compounds show similar properties and we 
will discuss ZnV204 mainly. At T st = 50 K, a cubic to 
tetragonal structural phase transition is observed |||9],|l0| , 
where the crystal is uniformly compressed along the c 
axis with c/a ~ 0.994 without any magnetic ordering 
HH . This lattice structure is just what we expect by 
taking 9 = 0, which leads to Q v > and Q u = in 
our theory. It is important to consider the effect of the 
difference of the u parameter (u = 0.260 |7|, |To| , |ll| ) from 
its ideal point (1/4), which exists even above T st . The 
distance between vanadium and oxygen, Zyo > is given by 
a(l — 4Z\m)/4, and cos/VOV = — 8Au in the lowest or- 
der of Au, where Au = u— 1/4. When the V-O-V angle 
is just the right angle, the superexchange path through 
the oxygen vanishes and the spin-spin interaction may be 
ferromagnetic (FM) due to the direct exchange. The ex- 
perimental results imply that cos _1 (— 8Au) is sufficient 
for the AF superexchange interaction to overcome the 
direct FM one. By the distortion due to the Q v mode, 
the cosines of the V-O-V in perpendicular and parallel to 
c axis change to — 8Au — Ap/ V6a and — 8Au + 2p/ \/6a, 
respectively. The change of the angles means that the 
coupling constant ^-\ Q is positive, see eq. (|l0|). There- 
fore for a positive Q v , as the experiments indicate, there 
is an enhancement of the AF exchange in the perpendic- 
ular direction to the c axis and \u) spin state is stabilized 
at each tetrahedron. 

The experimental magnetic susceptibility, Xj shows a 
weak temperature dependence above 100 K with a large 
Weiss constant, 8 W = 420 K fj. Below T clus = 95 K, 
X shows a splitting depending on the zero-field cooling 
(ZFC) and field-cooling (FC), which was interpreted as 
development of a short range or cluster ordering |^| . But 
its origin is not clear yet. Below T c i us , x shows only a 
gradual increase with decreasing temperature ||[7| and at 
T st , x drops with a sharp cusp structure This struc- 
ture may be consistent with the present scenario. Since 
the spin-singlet ground states themselves make no con- 
tribution to Xj the- spin part concerned with thermally 
activated triplet states and the Van-Vleck type orbital 
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part may be relevant to the total susceptibility. The 
separation of the doubly degenerate spin-singlet at each 
tetrahedron (denoted by As=o(T)) tends to suppress the 
spin part and the orbital fluctuation in the vertical plane 
to [111] axis is also suppressed. Both terms reduce \ with 
a sudden drop accompanied by the structural transition, 
which is qualitatively consistent with the experiments. 
With decreasing temperature, the increase of As=o(T) 
and the temperature dependent part of x, still slightly 
increasing around T st , may balance to result in the al- 
most temperature independent x below T st ||||] . 

At Tn = 40 K, a magnetic ordering sets in. The AF 
structure is observed by neutron diffraction experiment 
at 4.2 K (lO), where ID AF chains along [110] and [lTO] 
directions stack one after the other along the c-axis with 
the easy magnetization axis parallel to the c axis. This 
magnetic structure is consistent with the change of the 
AF coupling constants brought by the Jahn- Teller distor- 
tion of the Q v mode. To understand the AF long range 
order from the 3D VBS-like \u) spin states symmetrized 
on every vertex realized at T st , it is necessary to consider 
the effects of higher triplet states. 

Finally, let us discuss the electronic entropy, S e i, of 
ZnV2C>4. According to the trigonal distortion by the pos- 
itive Au, the threefold degenerate t2 g states are expected 
to split into the higher a\ g state and lower e g states. This 
scheme of the splitting is expected since the wave func- 
tion of the ciig and e g states are parallel and perpendic- 
ular to the [111] direction, respectively. On the other 
hand, for L1V2O4, the LDA+C7 calculation |TJ| suggests 
that a\ g state becomes lower than the e g to gain the 
coulomb energy of the singly occupied electron level of 
the (3d) 15 configuration, but this may not be the case for 
the (3<i) 2 state. Suppose that the orbital degree of free- 
dom is almost quenched at higher temperatures, then the 
main contribution to the electronic entropy around the 
structural transition is the spin degrees of freedom in the 
spin singlet sector, which is estimated to be R In 2 = 5.76 
J/mol-K and the lattice part additionally. This is consis- 
tent with the experimental value of S e i ~ 7 J/mol-K just 
above the T st |§. 

Another possible scenario for the structural transition 
is, of course, the usual Jahn- Teller distortion at T st . For 
this scenario, it is necessary to assume that the a\ g state 
is lower than the e g states. However, the entropy associ- 
ated with the orbital degrees of freedom is 2i?ln 2 = 11.52 
J/mol-K, which is much larger than the entropy observed 
experimentally. Note here that the spin degrees of free- 
dom are not included yet. It is also worth to mention 
that for the present scenario of the spin-driven lattice 
distortion, the energy scale should be a fraction of the 
exchange interaction J ~ 10 2 K. The smallness of the 
T st and structural metamorphosis at T st are compatible 
with this small energy scale. On the other hand, for the 
conventional Jahn- Teller effect concerned with the single 
electron states, the energy scale is typically a fraction of 
an electron volt (~ 10 4 K). 

In conclusion we have studied the frustrating quan- 



tum spin-1 system on the pyrochlore lattice as a theo- 
retical model of the insulating spin-1 spinel compound. 
In this problem, it is essential to lift the many-fold de- 
generate spin-singlet manifold to discuss the low energy 
physics. As such a mechanism, a magneto-elastic inter- 
action is considered in this letter and we find the sponta- 
neous breakdown of the lattice symmetry, which is con- 
sistent with the structural phase transition observed in 
ZnV204 and MgV204. Following this scenario, we can 
understand qualitatively the low temperature behaviors 
of the magnetic susceptibility and the electronic entropy, 
which seems to be difficult to be explained by the usual 
Jahn- Teller effect. It is natural to expect that spin-driven 
structural transition discussed in the present paper may 
be generalized to integer-spin pyrochlore systems with 
isotropic interactions. 
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